Optical absorption for parallel cylinder arrays 
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Abstract 

We study the long wavelength electromagnetic resonances of interacting cylin- 
der arrays. By using a normal modes expansion where the effects of geom- 
etry and material are separated, it is shown that two parallel cylinders with 
different radii have electromagnetic modes distributed symmetrically about 
depolarization factor ^. Both sets couple to longitudinal and transverse com- 
ponents of the external field, but amplitudes of symmetric depolarization 
factors become exchanged when considering longitudinal or transverse polar- 
ization. We also find that amplitudes satisfy sum rules that depend on the 
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ratio of the cylinders radii. 
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I. INTRODUCTION 



The optical properties of ordered cylinder arrays has become a subject of much recent 
interest mainly because of their potential use as photonic crystals. Jl]^] Theoretical re- 
search once done for spherical particles 0-0] has lately been applied to cylinders 



It is well known that in the long wavelength limit the optical properties of a dilute com- 
posite of microscopic spherical particles are well described by mean field theories such as 
Claussius-Mosotti or Maxwell-Garnett. These theories are essentially based on a dipolar 
approximation assuming that the particles are sufficiently far appart so that it is possible to 
neglect contributions from higher order multipoles. As particles become closer, however, this 
approach is no longer valid. Several models have been presented to overcome this difficulty, 
among wich the theory of normal modes has been shown to be particularly convenient since 
it makes possible to expand the response of the system in terms of optical resonance terms, 
where dielectric properties appear separate from geometrical factors. |lTT| , p!2| The simplest 
system exhibiting the effects of interactions is a pair of identical particles very close to each 
other. A pair of particles of the same material and form but different size is the simplest 
non-symmetric system of interacting particles 0. 

Recently, a model to study arbitrary cylinder arrays made of the same material has been 
constructed and applied in detail to a pair of identical cylinders. |]TU[ The response of periodic 
arrays of identical parallel cylinders including proximity effects is also easily treated within 
that formalism. The method follows a normal modes description appropriate to the long- 



wavelength limit first proposed by Bergmann [|lT],|12|, and makes use of a basis of cylindrical 
harmonics solutions to Laplace's equation. Modes are characterized by depolarization factors 
and strengths, defined in such a way that an isolated cylinder exhibits a depolarization 
factor I and unit strength. We here follow a similar procedure to study a pair of non- 
touching parallel cylinders of the same material but different radii, for different polarizations 
of the external field. We show that a difference in radii does not alter the property that 
depolarization factors are symmetric about ^, although in contrast with the equal radii case. 



3 



mixing of transverse and longitudinal modes occurs. This property is no longer valid when 
the dielectric function of the cylinders is different. We find that all normal modes are active 
for an external field perpendicular to the cylinders axis, whether parallel or perpendicular 
to the plane containing the axis. Modes inactive for identical cylinders have very small 
strengths when their radii are not very different. Strengths of modes with depolarization 
factors smaller than | are exchanged with those of depolarization modes bigger than |, when 
the direction of the electric field changes from parallel to perpendicular. 

In Sec. we get the multipolar moments and the absorption cross section for a pair 
of unequal cylinders. In Sec. |T| we present and discuss our numerical results. Finally in 
Sec. 1^ we summarize our conclusions. 



II. THEORY 

We consider a set of N parallel, infinite, uncharged cylinders of dielectric function 
placed in a homogeneous medium of dielectric function 62, excited by an external electric field 
whose wavelength is much longer than the cylinder radii or separation between cylinders. 
The charge distribution they acquire may be described in terms of individual multipole 
moments qmj obeying the equations, p 

Qmj = -amj {Vm+Yl ^mj Qm'j' ) , (l) 
m'j' 

Here m is a positive or negative integer labeling the angular momentum component along 
the cylinder axes, j = 1,2, is a particle index, amj = \m\ a^''"'(£:i — e2)/{si + 82) are 
the multipolar polarizabilities of cylinder cylinder j of radius aj, and Vm are the coefficients 
in the expansion of the external potential in terms of cylindrical harmonics. The coupling 
coefficients are given by, ||T0[ 

if m ■ m' > 

m m \m\+\m I ' 

I I I I pjj, 
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where {pjj', Ojj/) = pji — pj are polar coordinates in the x-y plane giving the relative position 
of cylinder j' with respect to cylinder j. 

As discussed in reference [|10|, if the cylinders are of the same material one can separate in 
Eq. (1) terms depending on the material susceptibility x from those involving the geometry 
of the array. We intend to follow here the same procedure and write Eq. (|I]) as 

J2ix~'5,,,+Hli')x,, =U, (3) 

where /i represents the pair of indices {niyj), and 

if™/ = 2n {6mm' Sjj> + |m mV^^ a!,'"' o}^ ^ A^j ) (4) 
fm, = -2n \m\'/' Vm, (5) 

^-^i ^ I u/') \m\ 

Note the important feature that matrix H depends on geometry only and its eigenvalues 
{4-7™^} define the depolarization factors {n^} of the array. For later convenience we write 
H = 27r(I + B), with I the unit matrix and -8^/ = 1"^ m'|^/^ a^"^' a}^ '^m/ ■> so that the 
depolarization factors {n^} and eigenvalues {A^} of B satisfy the relation, 

n, = \{l + K). (7) 

Because of the property -B^/ = if m ■ m' > (see Eq. (2)), we write rows and columns 
of matrix B with indexes m and m' following the sequence 1, 2, . . . , — 1, —2, . . ., resulting in 
matrix B written in terms of a new real matrix b of half its dimension, as 



B 



h 
b 



(8) 



From now on we use index m and m' as positive integers, and write the elements of matrix 
b as follows, 

C = - ^ , „ ' —, (9) 

""^ ^ ' m\m'\ 0^1' 
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It can be shown that the eigenvalues of matrix B come in pairs with opposite sign 
A^j = ±£^, where £^ are the eigenvalues of b (see the appendix |^). As follows from Eq. 
(7) the depolarization factors are then symmetric about the value 1/2. The components of 
vector Xfj_ can de written in terms of the eigenvalues of matrix b and elements of matrix u 
that diagonalizes b through 

u^^bu = i . (10) 

In the case of a uniform electric field Eq parallel to the plane containing the cylinder axes 
and perpendicular to the latter (parallel field geometry), Xmj = X-mj- Vector x+ = Xmj is 
then given by 

x+ = (u s-^u-^)f+, (11) 



where 



5w5jy(x" + 27r(l + i^j)), (12) 



fmj = SmlT^ttjEo. (13) 



In the case of an electric field perpendicular to the plane containing the cylinder axes (per- 
pendicular field geometry), Xmj = —X-mj- Vector x_ = x^mj is then given by, 

x_ = (ur-^u-^)r, (14) 

where now 

fmj' = 5mm'5jj'{.X'^ + 2vr(l - (.mj)), (15) 

fZmj = SmiiTrajEo. (16) 

For a pair of unequal parallel cylinders with radii ai and 02, and axis at a distance R we 
define dimensionless parameters (3 = a^jax and 5 = Rjax. Because of the properties Ix^j = 
and = (~/^)™~™'^m?5 we write rows (columns) of matrix b with particle index j (j') in 
the sequence 1, 2 resulting in matrix b written in terms of a smaller array g as follows, 



g 
g 

where g is the transpose of g, and the elements of matrix g are given by 



(17) 



9. 



I = {-r vW 



m\m\ 



(18) 



Resuhs for a pair of cyhnders with an external field in the parallel or perpendicular config- 
uration can be written in terms of a single normal modes expansion as 



E 



'--■mi J± 



(19) 



In this expression the upper (lower) sign corresponds to the parallel (perpendicular) con- 
figuration, with (m',j') labelling the excitation modes of the pair as a coupled system. It 
gives just half of the multipoles; the others are obtained from the corresponding symme- 
try property as given in the paragraph preceding Eqs. ( pT]) and (p!!]). We have defined the 
depolarization factors of modes, 

1 



n± 



m' f 



-1^1 ± im'j') , 



and coefficients corresponding to strength of modes. 



(20) 



C, 



(21) 



We have also defined 



/+ — /iti 

/- = /ri,i . 

We find that coefficients C^/ satisfy the sum rules 



^ml — "ml 

m'j' 

E ^rn2 = P^ml 



(22) 
(23) 



(24) 
(25) 
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It can be shown that, as with the original matrix B, eigenvalues {imj} of matrix b come 
also in pairs with opposite sign; therefore the sets of depolarization factors {n+mj} and 
{n^mj} are identical (see the appendix). A given depolarization factor exhibits a different 
strength depending on the direction of the external field. As seen in the normal modes 
expansion given by Eq. (^) the same strength coefficients C^^f appear for depolarization 
factor n+rn'j' in the parallel field response and for n_m'j' in the perpendicular field response. 
Then, strenghts corresponding to depolarization factors symmetric around value 1/2 are 
exchanged between responses corresponding to fields parallel or perpendicular. 
The magnitude of the electric dipole moment for the pair can be written as 

2 \ " '-^11 "T P^12 771 fna\ 



m'j' ■ -tm j 

where p+ corresponds to a parallel (perpendicular) field. The absorption cross section is 
proportional to the imaginary part of the factor accompanying Eq in the previous expression, 
a quantity we identify as the complex effective polarizability of the pair. Thus we arrive at 
a normal modes decomposition for the absorption cross section of two parallel cylinders, 

a± ~ Imaginary{ '\ '\, } (27) 

m'j' A ' lt'±m'j' 

where (T+ (cr_) corresponds to a parallel (perpendicular) field. Notice that, as follows from 
Eq. (pID, the numerator in the previous sum can be written as 

K?' + /9</)', (28) 

and is always positive definite. According to the sum rules given by Eqs. (^) and (^), the 
sum of the numerators in expansions and (^) is 1 + a feature we use in calculating 
the normalized strength of modes. 



III. NUMERICAL RESULTS 



We have solved numerically the eigenvalue equation for matrix b for the case of a pair of 
parallel cylinders of radii ai and 02, and have calculated the depolarization factors n^'j' and 
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strength coefficients according to Eqs. (pO]) and (pTl). We have studied in detail the 

normal modes expansion for the dipole moment of the pair as given by Eq. (p6D. Our most 
important finding is that when the radii are not equal, modes with depolarization factors 
above and below the value | mix for all orientations of the external field. This is known not 
to happen when cylinders are equal. flO 



In obtaining numerical results we use the dimensionless parameters a = R/{ai + 02), 
that measures the center to center distance, j3 = ci2/fli, measuring how disimilar the radii 
are, and fi = {R — ai — a2)/o.ii measuring the border to border distance. In Fig. |I] we plot 
modes for very close cylinders (a = 1.10) with one radius three times bigger than the other 
{(3 = 3). Modes are for the parallel configuration, while those for the perpendicular case 
are obtained by mirror reflection about depolarization factor |. Note that the placement 
of modes are symmetric about this central value so that, as far as position is concerned, 
they are indistinguishable in both configurations. The figure shows the modes with largest 
amplitude, while weaker modes accumulate around n = ^. The sum of amplitudes is 10, as 
required by the sum rule mentioned after Eq. (pSf). Labels, included in order to match with 
labelling in Fig. ^ are arbitrary. 

Figure H shows depolarization factors (a) and normalized strengths (b) in terms of the 
parameter j3, at fixed a = 1.1. The sum over all strengths is unity, having been normalized 
by the factor l+/3^. Thus the results at /? = 3 correspond to amplitudes shown in Fig. ^with 
the normalization factor 10. In changing j3 we keep constant the parameter a by changing 
the center to center distance R accordingly. Note that modes with essentially zero strength 
at /5 = 1 become important when increasing this ratio. Labels are arbitrary and are used 
just to relate the data in different figures. 

Results shown in Fig. ^ were obtained by changing (3 and the center to center separation 
R, but keeping constant the border to border distance at the fixed value /i = 0.4. In this case 
we note that the depolarization factors move appart with increasing (3, while at constant 
a (Fig. 2 (a)) they get closer. This is because the relation between the edge to edge and 
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center to center parameters is /i — {a — + P), indicating that as /3 goes to infinity, so 
does if (J is kept constant. Thus all modes should converge to the isolated cyhnder value 
n = 1/2 in this case, while if is kept constant, the modes converge to those of a cylinder 
in front of a plane at the same distance. 

IV. CONCLUSIONS 

In summary, we have shown that the absorption cross section of a pair of parallel cylin- 
ders of the same material but different radii contains modes whose depolarization factors 
are symmetrically distributed around | with amplitudes depending on the direction of the 
external field. When the field changes from the parallel to the perpendicular configuration, 
amplitudes corresponding to symmetric depolarization factors about the central value 1/2 
are exchanged. 
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APPENDIX A: 



We consider the eigenvalue equation for operator T, 

Tx = Ax , 

in a basis of dimension 2M. Here T has the form, 



(Al) 



A 
B 



(A2) 



where A, B each has dimension M. In writing the eigenvectors x in terms of two smaUer 
vectors w and v of dimension M the eigenvalue equation is cast into the form, 



(A3) 



or 






A 




w 


= A 


w 


B 







V 




V 



Av — Xw 
Bw = Xv . 

Prom there we get the separate eigenvalue problems, 

{AB)w = X'^w 
{BA)v = A^^; , 



(A4) 
(A5) 



(A6) 
(A7) 



both having the same eigenvalues A^. In solving for the corresponding eigenvectors and 
writing them as columns we get matrices w and v that diagonahze matrices AB and BA. 
They can be used to form a matrix U as. 



U= ^ 



w w 

V —V 



(A8) 



which diagonalizes matrix T according to the relation 

U-iTU = A , 
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(A9) 



with matrix A given by 



(AlO) 

Here A (—A) is a diagonal matrix formed by the positive (negative) square root of the 
eigenvalues of matrices AB or BA. Therefore the eigenvalues of matrix T come in pairs 
with opposite sign. In the case of a pair of cylinders, we find the previous feature two times. 
It first happens because coupling coefficients A!^^ are zero if m and m' have the same sign, 
and then occurs also because they are zero for j — j'. The dimensionahty (4M)X(4M) of 
the original eigenvalue problem is seen to be reduced to MXM dimensions. 




12 



REFERENCES 

[1] E. Yablonovitch, T. Gmitter, K. Leung, Phys. Rev. Lett. 67, 2295 (1991). 
[2] N. A. Nicorovici and R. C. McPhedran, Phys. Rev. E 54, 1945 (1996) 
[3] K. Ohtaka, T. Ueta, K. Amemiya, Phys. Rev. B 57, 2550 (1998) 

[4] W. T. Perrins, D. R. MacKenzie, and R. C. McPhedran, Proc. R. Soc. London A 369, 
207 (1979) 

[5] F. Clare, Phys. Rev. B 30, 4989 (1984); 35, 406 (1987) 

[6] R. Rojas and F. Claro, Phys. Rev. B 34, 3730 (1986) 

[7] I. Ohvares, R. Rojas and F. Claro, Phys. Rev. B 35, 2453 (1987) 

[8] A.V. Radchik, G.B. Smith, A.J.Reuben, Phys. Rev. B 46, 6115 (1992). 

[9] A.J. Reuben and G.B. Smith Phys. Rev. E 58, HOI (1998) 

[10] R. Rojas, F. Claro and C. R. Proetto, Phys. Rev. E 62, 5688 (2000) 

[11] D.J. Bergman, Phys. Lett. 43C, 377 (1978) 

[12] D.J. Bergman, Phys. Rev. B 19, 2359 (1979) 

[13] Lord Rayleigh, Philos. Mag. 34, 481 (1892) 

[14] R.C. McPhedran, D.R. McKenzie, Appl. Phys. 23, 223 (1980) 

[15] R.C. McPhedran, W.T. Perrins, Appl. Phys. 24, 311 (1981) 

[16] J.B. Keller, J. Math. Phys. 5, 548 (1964). 

[17] J.M. Pitarke, F.J. Garcia- Vidal, J.B. Pendry Phys. Rev. B 57, 15261, (1998). 



13 



FIGURES 

FIG. 1. Mode amplitudes for a pair of unequal cylinders under a uniform electric field in the 

parallel field configuration. The raddi ratio equals /3 = 3, and the separation parameter is a = 1.1. 

FIG. 2. Depolarization factors (a) and normalized strengths (b) as a function of size parameter 
/?, for a pair of unequal cylinders with fixed separation parameter cr = 1.1 under a uniform electric 
field. Labels correspond to those in Figure 1. 



FIG. 3. Same as Fig. 2, but with fixed border to border parameter /x = 0.4. 
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